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TIME -DEPENDENT DOWNWASH AT THE TAIL AND TEE PITCHING 
MOMENT DUE TO NORMAL ACCELERATION AT 
SUPERSONIC SPEEDS 
By Herbert S. Ribner 

SUMMARY 


The time -dependent downwash behind a wing in a supersonic stream 
is analyzed for the case when the angle of attack varies linearly with 
time. The result is applied to the calculation of the contribution of 
the horizontal tail to the pitching moment and lift due to normal accel- 
eration of the airplane. The method employs an extension of an unpub- 
lished solution of the linearized potential "equation for unsteady flow 
by Clifford S. Gardner. The pitching moment due to normal acceleration, 
together with the damping in pitch, determines the damping of the short- 
period mode of longitudinal oscillation for an airplane. 


INTRODUCTION 


The investigations of Garrick and Rubinow (reference 1) and others 
have shown that a two-dimensional wing may experience certain unstable 
torsional oscillations at low supersonic speeds. These oscillations are 
in pitch without coupled vertical oscillations of the center of gravity. 
The same behavior may also occur with three- dime nsional wings. More 
generally, coupled oscillations in these two degrees of freedom, termed 
"short -period" oscillations, may be expected. The damping of these 
oscillations is governed by the sum of the damping-in-pitch deriva- 
tive Cm^ and the pitching-moment coefficient due to the normal-acceler- 
ation derivative C^, where a is the time rate of change of angle of 
attack. The sum is therefore of great significance, inas- 

much as its sign determines whether the motion is stable or unstable. 

An unstable sign usually arises from the component C^. 

The derivative for an airplane is compounded of a contribution 

from the wing and another from the horizontal tail. At low subsonic 
speeds the wing contribution is a negligible factor, but at supersonic 
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speeds it may "become significant. The relevant theory for the wing at 
supersonic speeds and some calculations for specific plan forms are 
given in references 1 to 7. 

The contribution of the horizontal tail to Cm^ depends to a large 

extent on the time -dependent downwash field of the wing when the angle 
of attack varies with time. Standard techniques based on the solution 
of Laplace's equation (that is, the assumption of incompressible flow) 
have been used for low subsonic speeds. These methods fail at higher 
speeds when the compressibility of the air must be taken into account. 
Then this time -dependent downwash, in the theory of small disturbances, 
must satisfy the time -dependent form of the- Prandtl-Glauert equation. 
(See equation (l) herein.) In the present paper a solution of this 
problem is obtained for. supersonic speeds, and some considerations are 
given for subsonic speeds. 

The evaluation of the time -dependent downwash employs an exten- 
sion of an unpublished paper by Clifford S . Gardner of New York 
University. (The relevant part of his work is set forth in an erratum 
sheet that was issued for reference 6 and is repeated in the body of 
the present paper.) It. is shown that the problem of calculating the 
pressure distribution due to the time rate of change of angle of 
attack a can be reduced to the well-known problems of calculating 
the pressure distributions due to steady angle of attack and to steady 
pitching. Gardner's velocity potential applies only ahead of regions 
affected by the trailing edge of the wing. The necessary modifica- 
tion for these regions, particularly the region behind the trailing 
edge, is considered herein. 

The time -dependent downwash obtained in this manner is applied to 
the calculation of the contribution of the horizontal tail to the 
normal -acceleration derivatives Cm& and Cl^. The derivation is 

limited to the case where the tail lies in the chord plane of the wing 
(z = 0). In addition, formulas for the contribution of the wing to 
these derivatives, evaluated from Gardner's work, Eire given for 
convenience . 


SYMBOLS 


x, y, z rectangular coordinate system fixed in horizontal wing: 

x-axis chordwisej y-axis spanwisej z-axis upward 


t 


time 


V 


stream velocity, directed parallel to x 
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Cm 


speed of sound in free stream, 
stream Mach, number (V/a) 

■velocity potential 
angle of attack 

(assumed constant herein) 

angular velocity of pitching 
angular frequency (sinusoidal motion) 
value of 0 for q = 1, a = a = 0 
value of 0 for a = 1, q = a = 0 

part of X associated with trailing vortex system of wing 

X-,, y~, and z -velocity components of specified potential 
flow (see text) 

upwash due to angle of attack (aX z ) for z = 0 

upwash due to pitching (qt z ) for z = 0 

local pressure on lower surface of airfoil minus local 
pressure on upper surface (Lift /Unit area) 

density of air 

lift 

pitching moment 

/ T. \ 

lift coefficient 
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pitching -moment coefficient 
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airfoil area 

airfoil mean aerodynamic chord 
airfoil incidence relative to x-axis 
x-coordinate of reference point in tail (tail arm) 
local sweep angle along trailing edge (A(y)) 


Fl(x,y)' = 1 + ^ 


“o 

? 2 (x,y) ■ — g 

Fl(2) 

f 2 (z) 

w(2) 

V !) 


KS) - x (&) 


x te\S 


2 cos 2 A - 1 


Mq 2 ^ IE aV 


dx 


VMq cos A 

average of F-^x^y) over tail span for x = Z 
average of F 2 (x,y) over tail span for x = Z 
average of w(x,y) over tail span for x = Z 
average of w^(x,y) over tail span for x = Z 
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Subscripts : 

TE 

T 


wing trailing edge 
tail (horizontal) 


When x, y, z, or t are used as subscripts, the respective 
partial derivative is indicated. For example, 

<k 

™ ‘ 3k 

^ " a* at 


ANALYSIS 


The linearized partial differential equation for unsteady- 
supersonic flow is 


- f*yy - + % #rt + Jg At - 0 (1) 

The boundary condition on a wing experiencing a constant time rate of 
change of angle of attack is 

0 Z = -aVt (z — >0) (2) 

The disturbance potential for snail disturbances must satisfy both 
equations (l) and (2). 

In an unpublished paper Clifford S. Gardner has shown, in effect, 
that a suitable solution is 


0 MO 2 . 

a b 2 


t - 


Mo 2 xN 
VB 2 J 


X 


( 3 ) 


where + is the steady-state potential corresponding to a unit 
pitching velocity about the y-axis, and X is the steady-state 
potential corresponding to a unit angle of attack. (See also the sec- 
tion entitled "Discussion.") The fact that equation (3) is a solution 
can be verified by direct substitution into equations (l) and (2). 
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Thus, Gardner has shown that the time -dependent potential for an angle 
of attack at may be compounded of two time -free or steady- state 
potentials, one for a constant angle of attack, and the other for 
steady pitching. (Equation (3) is limited in direct application to 
flow regions uninfluenced by the trailing edge. This limitation will 
be elaborated later.) The first-order treatment of sinusoidal oscilla- 
tions in reference 5 (equation (19) and the interpretation therein) 
leads to a relation corresponding to equation' (3) . The relation in 
reference 5 is derived, however, as an approximation for slow oscilla- 
tions; whereas Gardner's relation is exact for a constant accelera- 
tion (a = constant). Equation (3) is also contained, in effect, as a 
special case of a solution for a more general unsteady motion given in 
reference 8. The plan forms for which the derivation in reference 8 
is applicable are not completely general in contradistinction to those 
considered herein. 


The lift distribution at time t = 0 for the angle of attack at 
is obtained from the upper-surface potential by the linearized 
Bernoulli equation for unsteady motion, 

AP = 2p(V0 x + <h) z=+0 (4) 


AP = 


AP = 


2pVife- t x 


. JL) 

VB2 * VB2 / 


Mo 2 x, 


^(MO (AP)^! - 2 P X 


(5) 


where (AP) _ 1 is the lift distribution for unit steady pitching 
velocity about y-axis, (AP)^^ is the lift distribution for unit 

steady angle of attack, and equation (4) contains the implicit 
assumption that (0 ) z __q = ~(P) z =+0’ ® ie c ^°^ ce t = 0 elimi- 

nates the lift due to angle of attack and leaves only the increment 
due to time rate of change of angle of attack. 

Integration to obtain the lif-trand moment and reduction to 
coefficient form yields 
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( 6 ) 


c m& = 





18 + d eJMfh 


(7) 


where the integrations are carried over the wing plan form. 

Equations (6) and (7) provide an evaluation of the acceleration 
derivatives Cl^ and for the wing in terms of quantities that 

are presumed to he known for the pitching wing and the wing at an angle 
of attack. The application of these equations is limited to wings with 
supersonic trailin g edges. (The component stream velocity normal to 
the trailing edge is supersonic.) The reason for this limitation is 
shown later. 


Time -Dependent TJpwash at Tail 

Wings with supersonic trailing edges. - For the evaluation of the 
contribution of a horizontal tail to and it is first 

necessary to evaluate the instantaneous upwash velocity 0 Z at the 
tail location. Differentiation of equation (3) would afford this 
quantity if equation (3) were still valid behind the wing trailing 
edge. Equation (3) as applied behind the trailing edge still satisfies 
the differentia l equation (l) and the boundary-condition equation (2) 
(limited to the wing surface) if the values of ^ and X therein 
assume values appropriate to the region behind the wing, including the 
contribution from the trailing vortex sheet. Equation (5) , however, 
shows that the lift fails to fall to zero there as it should: the 
terms in $ x and X^ with i and X as previously defined, are 
known to vanish but the remaining term in x does not vanish. 

Thus equation (3), in its present form, leads to a spurious lift 
behind the wing of amo unt -2pXa/B^ per unit area. The situation may 
be corrected, however, by an application of Lagerstrom's concept of 
cancellation of lift. There is superposed behind the wing an additional 
flow Uj_, Vj_, w-j_ that gives rise there to an equal lift of opposite 

sign. Note that this lift is independent of time. The value of uq 
along the upper surface (z = +0) of the trailing vortex sheet is then 
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uq(x,y,+0) = ^ (8a) 

according to the steady-state linearized Bernoulli equation, which is 
equation (4) with 0^ = 0. The further evaluation of uq and W]_ 

(vq is not needed) is limited at first to the relatively simple case 
where the trailing edge of the wing is supersonic. 

Equation (8a) and its counterpart for the lower surface constitute 
a boundary condition on uq. The distribution of uq in space (z / 0) 
must be such that the corresponding potential satisfies the steady- 
state linearized potential equation, which is equation (l) without the 
time -dependency terms. (This potential, being time free, will also 
satisfy the complete equation (l).) By differentiation of equation (l) 

50 

with respect to x, U 1 = r— is seen to satisfy the Bame partial 

differential equation as 0 does. The part of X,Xq r that can be 
ascribed to the trailing vortex system of the wing satisfies the same 
boundary condition at z = ±0 (with the factor a/VB^) and the same 
partial differential equation. The solution can be shown to be unique, 
and thus uq may be identified with Xtyd/VB 2 everywhere without 
restriction to z = ±0. Behind the wing X = X^ r for z - ±0 but 
does not for other values of z. Equation (8a) is thus generalized to 

uq(x,y,z) - (8b) 

In what follows, however, the small difference between X^ and X 
behind the wing is neglected to simplify the calculation. 

Associated with this additional u-velocity uq is an upwash 
velocity wq. This upwash may be determined from the irr ot at io n a l ity 
condition 

dwq 8uq 

ST = 3T 

By virtue of equation (8b) with X in place of %fc r this upwash is 

^1 a y 
= ;? * 

For z = 0 there results upon integration 

n = + 4? * ** 


(9)* 

( 10 ) 4 
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The quantity wq^qg is the upwash induced just behind the trailing 
edge of the wing (leading edge of the superposed flow) by the superposed 
flow. This edge has been assumed supersonic; hence, in its vicinity 
the Ackeret two-dimensional flow equations may be applied to the flow 
components normal to the edge (that is, simple sweep theory may be 
applied) . The following equation is obtained: 

R.P.\/mo^ cos 2 A - 1 

W 1,TE = -U 1,TE3 cosT ^ Z = + 0 ) 

where A = A(y) is the local angle of sweep of the trailing edge. 

(The real part designation R.P. is superfluous here but has a later 
use.) By use of equation (8a) this relation may be incorporated in 
equation (10) in the form 


w n = 


, <X / R.P.n/Mo^ cos 2 A - 1 r x v 

w 1 (x,y,0) = — r -X™ — + / \ 

VB 2 V ^ COS A J rjTg 


dx 


(H) 


z=+0 


Equation (ll) gives the upwash component wq of the 
flow uq, vq, wq that was superposed to cancel the spurious lift 

behind the wing. This equation applies to points behind the wing and 
in the same plane. The upwash wq is to be added to that computed 
as j z from equation (3). 

Wings with subsonic trailing edges . - The foregoing development 
for wq has been limited to wings with supersonic trailing edges. A 
rigorous determination of the superposition flow uq, vq, wq for wings 
with subsonic trailing edges is much more difficult and the process is 
only indicated herein. For the wings with subsonic trailing edges, the 
steady potential flow uq, Vq, Wq to be superposed must satisfy the 

boundary condition, equation (8a), (and its counterpart for z = -0) 
behind the trailing edge, and must, in addition, satisfy the condi- 
tion wq = 0 on the surface (z = ±0) ahead of the trailing edge within 
the region of the trailing -edge disturbance. A flow of this kind may 
itself be built up by superposition of simpler "mixed-wing" flows in a 
variant of the manner elaborated in reference 9. 

A rigorous evaluation of wq may not be necessary for practical 
applications, however, if the region of the wing influenced by subsonic 
trailing edges is not large. A suitable approximate value of wq may 
be obtained by a modification of the procedure leading to equations (10) 
and (11 ). Equation (8a) relating the surface value of uq to the 
surface value of x is still true for wings with subsonic trailing edges. 
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The generalization for values off the surface (equation (8b)), however, 
no longer holds rigorously. The assumption is now made that an equa- 
tion of the form of equation (8b) holds approximately with X in place 
of X tr . Then it follows that, to this approximation, equation (10) 
applies to wings with subsonic trailing edges. 

For this application, the value of wq at the trailing 
edge, Wq rjTg, is taken to be zero. This fact follows from the boundary 

condition wq = 0 ahead of the wing trailing edge, together with the 
observed fact that— wq must be continuous across a boundary on one 
side of -which- wq is specified and on the other side of which uq is 
specified, both of these velocities being finite. Equation (11) may 
now be recognized to apply for both subsonic and supersonic trailing 
edges because of the restriction of - the radical to its real part (R.P.). 
Thus, the term containing the radical, which is w q_ qzE, automatically 

vanishes for subsonic trailing edges. 

Remarks on wings in a completely subsonic stream . - The partial 
differential equation (1) applies equally well at subsonic and super- 
sonic speeds. Thus, equation (3) as applied to subsonic speeds is 
still a solution of equation (l) that satisfies the boundary condition 
of equation (2) on the wing but yields a spurious lift in the wake. 

Again an additional flow Uq, Vq, wq is required to cancel this spurious 
lift- behind the wing. The boundary conditions for this cancellation 
flow are specified as for wings in a supersonic stream with subsonic 
trailing edges, but - in the present case the region of trailing-edge 
disturbance covers the entire wing. Accordingly, the incremental lift 
corresponding to uq affects the entire wing surface. Equations (5) 
to (7) do not include this important contribution. 

The values of 'if and X specified in equation (3) are potentials 
obtained from lifting- surface theory. Expressions obtained from 
lifting-line theory, however, may be applied in the region behind the 
wing. These expressions should serve for the calculation of the time- 
dependent up wash at the tail location. 

Angle of attack at the tail .- The angle of the local flow relative 
to the tail chord is compounded of the tail incidence irp, the airplane 

0 Z wq 

geometric angle of attack &t, — and the upwash angle — + 

.. 0z W 1 

Ct<J! = lrjl + at + -y" + -y“ 

The component 0 Z is obtained by differentiating equation (3) (any z). 
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and, for z = 0, wq is given by equation (ll). In the following 
discussion, consideration is limited to the case where the tail lies 
in the chord plane of the wing (z = 0). Then, using equations (3) 
and (11), 


• / ^z' 

cifji = i^ + at 1 + 


'z=0 





^ cos^ A - 1 


VM cos A 





( 12 ) 


The term Xg in equation (12) may be identified as the upwash due 
to a unit wing angle of attack, and the term may be identified as 
the upwash due to a unit wing pitching velocity. It will be convenient 
to replace X z and \jr z by means of the definitions: 



(13) 


Then 




XqbR.P. x/Mo 2 co S 2 A - 1 

— + 

VMq cos A 


(14) 


Abbreviate this to 


a<j = irp + dt F 1 (x,y) + a, F 2 (x,y) 


( 15 ) 


Tail Contribution to G nd 

Equation (l4) or equation (15) gives the instantaneous angle of 
attack at any point (x,y) of a tail located in the z = 0 plane 



12 


MCA IN 2042 


relative to a system of axes with origin in the wing. The effect of 
the variation of aq> in the y-direction may he presumed for the 
present purpose to he adequately taken care of hy forming the arithmetic 
average over the span of the tail. Denote this average hy means of a 
har above Fq and T?q and above w and w^ in F^_ and F 2 . The 
variation of aq in the x-direction is more significant; - For practical 
application, however, it will suffice to use, iru effect, an average in 
the x-direction also. Thus, the angle of attack 

cc^i = iij + at F]_(2) + a F 2 (2) (z = 0) (l6) 

where x is taken equal to the tail arm 2, is presumed to apply 
uniformly over the entire tail area. The tail arm 2 is measured 
from the origin to same reference point on the tail, such as the tail 
center of gravity or the elevator hinge line. 

Equation (l6) foims the basis for obtaining the first approximation 
to the tail contribution to the derivatives and C^. Thus, the 

tail lift at, time t = 0 may be written 


^ = if (<^)t=o + ^(“rjt-o 

ScCm 

where dr = ^jr~« The choice of time t - 0 eliminates the lift due 

to the airplane angle- of attack and leaves just the increment due to 
time rate of change of angle of attack. By virtue of equation (l6) 
the tail lift may be written 



Then, upon differentiation with respect to d. 


dLq dLq — dLq _ 

sr = 555; Fl + 3^ F2 


(17) 






w. 


The pitching moment of the tail about its own reference axis is assumed 
in the present approximation to be small compared with the wing pitching 
moment. The only moment derivative with respect to ct contributed by 
the tail is therefore the component 


A 


3M = _ z 

3d " da 


1 
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Also, the & lift derivative contributed is the c omp onent 


. 5L &-T 

A S& = SST 

In nondimens ional coefficient form, these several equations become 

% ‘ ( CL a) T + +'<**) 


tt ■ («4 1 2 + K) T ^ < i8 > 

2VStZ ^CLrn 

^2 5dT ( 1 9) 

2VS T SCLrp 

AC ^ = ~sT “8T ( 2 °) 


Note that (ci^) t and (Cj-^ are merely the values of CLq. and Ql^ 

computed for the isolated tail considered as a wing in an undisturbed 
flow. The lift-curve slope Ci^ may be computed or estimated by well- 
known methods, and a formula for the evaluation of C]^ is given in 
equation (6). The functions F]_ and F2 are defined in equations (13) 
to (l6). The upwash parameters w/aV and w^yqc therein may be 

plotted once and for all for a given wing for use in calculations of 
the present type. 


DISCUSSION 
Downwash Charts 


The contribution of- the horizontal tail to the derivative Cj^ 
is seen to depend on a knowledge of the wing downwash due to angle of 
attack. -w/aV and of the wing downwash due to pitching -Vqjqc. 
(Furthermore , the tail contributions to the better -known deriva- 
tives Ch^ and Cm^ depend, respectively, on the same two quantities.) 
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Calculations and charts of -v/dV are already available for several 
types of wings at supersonic speeds (references 10 to 15 ), but data 
on w qjqc or its average over the tail span Wq/qc are lacking. This 
quantity is small compared with w/aV and the following simple con- 
sideration is sufficient for its estimation. 

It is easy to show that far behind the wing (compare references 10 
and 11) 


v q _ w 

qc ~ 2C La aV 


(x = °4 


( 21 ) 


if the shape of the span loading curves for pitching and angle of attack 
are the same. For differing shapes, equation (21) will at leastr provide 
the order of magnitude of Vqyqc. The ratio Cl ^20]^ is usually only 

a fraction of unity: for rectangular wings pitching about the midchord 
line the maximum value is l/6j for triangular wings pitching about- 


the --chord line a representative value is l/8. 


In the expression 


for Fg(z) (equations (l4) and (15), with bars added) w^yqc is 

multiplied by c and w/aV by x = l. Since the tail arm l is 
ordinarily two or three times the chord cT, the term in Wq/qc is seen 
to be quite small in comparison with the term in w/aV. Accordingly, 
the term in w/qc may frequently be neglected. In any event, equa- 
tion (21) should provide a sufficiently good approximation. 


Some Limitations on the Results 

The final equations (l4) to (20) are not rigorous. The approxima- 
tion was introduced in the evaluation of the upwash wq associated 
with the cancellation flow for the spurious lift - in the wake. The 
approximation lies in the integral that occurs in Fq as defined in 
equation (14). The simulation of the exact integrand by w/aV was 
made there to render the calculations more tractable. For wings with 

8 / x tr\ 

supersonic trailing edges, the exact integrand is » this 

integrand is believed to differ only slightly from v/aV. For wings 
with subsonic trailing edges, the exact-integrand has not been obtained 
because of the difficulty. The exact integrand is believe cLnot to 
differ greatly from w/aV, provided the trailing-edge disturbance does 
not envelope a large proportion of the wing surface. 

An error also occurs in equations (6) and (7) for Cq^, and 
for the wing when the trailing edges are subsonic. The magnitude of 
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the error likewise depends on the proportion of the wing surface 
enveloped by the trailing-edge disturbance. 

The method of derivation used herein employs the assumption that 
the angle of attack varies linearly with time. The airplane, however, 
may be expected to experience sinusoidal variations of angle of attack 
with time, and a sinusoidal variation is assumed in solving the equa- 
tions for dynamic stability. The values of C ®<x and Cl^ obtained 

herein are believed to be satisfactory approximations for sinusoidal 
motion so long as the square of the reduced frequency cuc/2V is small 

compared with B^/M^. This result is inferred from the fact that 
equation (3) is still approximately true for slow sinusoidal motions. 
This -result was proved in reference 5 (equation (19) and following 
discussions therein) . An examination of the terms omitted therein in 
the power expansion in cjd leads to the criterion for co previously 
noted. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., December 2, 194-9 
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